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We propose an improved viscosity model accounting for experiments of emulsions of two 
immiscible liquids at arbitrary volume fractions and low shear rates. The model is based on 
a recursive-differential method formulated in terms of the appropriate scaling variable which 
emerges from an analysis of excluded volume effects in the system. This variable, called 
the effective filling fraction, incorporates the geometrical information of the system which 
determines the maximum packing and reduces to the bare filling fraction for infinitely diluted 
emulsions. The agreement of our model for the viscosity with experiments is remarkable for 
all the range of volume fractions and viscosity ratio. 

I. INTRODUCTION 

Due to the central role that they play in many technological processes, the rheology of so^ 



.ogyof solid- 



liquid suspensions is a subject for which a large amount of work has been produced 
^, 7|. However, the rheological properties of emulsions of immiscible liquids have received much 
ess attention, despite the fact that they are also very important in many industrial applications 



3s attention, ciesp: 
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Emulsions present an interesting rheological behavior with characteristics similar to those of the 
solid-liquid suspensions. In particular, the relative viscosity rjr {(j)) of emulsions of nearly spherical 
droplets also diverges at certain critical value of the filling fraction: -q^ {4> i'c) oo. Also similar 
is the fact that for very dilute emulsions where the interaction between neighboring drops is absent, 
the relative viscosity follows the Einstein's like relation 

where the viscosity ratio K = ija/rjc contains the viscosity rja of the dispersed phase and the 



viscosity ijc of the continuous phase. Equation ([T]) has been originally introduced by Taylor 
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and is called Taylor's equation. The difference between Einstein's and Taylor's expressions lies 
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in the fact that in the later relation the coefficient multiplying the filling fraction (j) incorporates 
information about the nature of the dispersed phase through its viscosity r/^. However, when the 
viscosity of the dispersed phase is much larger than the viscosity of the continuous phase, X — s- oo, 
one recovers Einstein's expression r/ {(f)) = r^c (1 + 2.50). 

Although of fundamental importance, Taylor's equation does not reproduces the behavior of 
concentrated emulsions since in this case interactions between droplets becomes very important. 
Several efforts have been done in order to incorporate these interactions 1^, [3]. Among them we 
can mention the viscosity equation developed in Ref. based on a differential effective medium 
procedure 

(2) 

that takes into account the effect of the viscosity ratio K. A shortcoming of this expression is that 
it fails to describe the experimental data adequately at large concentrations. Pal improved this 
model by incorporating "crowding effects" through a critical filling fraction (j)c, sl: 



2 + hK J 

This model (called model 2) considerable improves the quantification of the rheological properties 
of emulsions of two immiscible liquids even at high filling fractions. However, comparison with 
experiments shows that it underestimates the value of the viscosity at intermediate volume fractions 
Sj. Pal also proposed the following model giving excellent agreement with experimental data at 
large volume fractions [q] 

however, it does not reduce to the correct Taylor's expression at low concentrations. 

Having in mind all these considerations, in this article we propose an improvement to the 
differential effective medium model that gives the correct divergence at the critical concentration, 
produces a better description of experiments and reduces to the right Taylor's expression at low 
concentrations. 



II. IMPROVED DIFFERENTIAL VISCOSITY MODEL: SCALING AND EXCLUDED 

VOLUME EFFECTS 

In its original article, Taylor calculated the viscosity of a dilute suspension of fluid drops under 
three main assumptions: i) The drops are small enough to maintain a spherical shape due to 
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surface tension, ii) no slipping exists at the interface between the drops and the host fluid and in) 
tangential stresses are continuous at the surface of the drop. Then, using hydrodynamic results 
for the velocity field out and inside of the sphere, he obtained the corresponding expressions for 
the components of the stress acting across unit area of the spherical surface and imposed on 
the tangential components the continuity condition. The normal stress is not continuous at the 
surface, fact which is ultimately related to the surface tension. As a result of these impositions, 
Taylor obtained the explicit expressions for the four constants appearing in the formulas giving the 
components of the stress tensor. Then, following Einstein's arguments, Taylor showed that in the 
case of liquid droplets, the factor 5/2 that multiplies the volume fraction (p in Einstein's relation 
Vr (0) = [1 + (5/2) (j)] must be replaced hy F = —5/2 [ijii + {2/b)rjc\ / [rjd + r]c\. This replacement 
leads to Eq. ([I]). 

The procedure followed by Taylor is very interesting since it emphasizes two different ingredients 
entering into the suspended phase correction Ar/ (= r]r — 1) in Eq. ([T|) 

The first ingredient is that the boundary conditions at the surface of the drop determine the 
stresses applied on it and thus control the contribution due of a single drop to the total stress 
of the composite system: particle plus continuous phase. The second ingredient comes from the 



consideration of the contribution of N droplets to the composite system. Following Landau 
one may argue that the total contribution by the drops 11^ to the total stress tensor comes from 
an average over the volume V of the system in the form 

nl = f / nW<iF. (6) 

where we represented the single drop contribution to the stress tensor by n^^"* which is proportional 
to F and with a the radius of the drop. The upper V in Eq.([6]) stands for the volume average 
and the factor N accounts for the contribution of the independent particles. However, this 
average is strictly valid only when the system is made of point particles. This restricts the validity 
of Eq. d^D to very diluted systems. 

Thus, if one considers that a drop has a volume then the average must be performed over 
the free volume accessible to the particles, which is defined by: Vjree = V — cNVd- Here, Vd is 
the volume of one particle and c is a constant taking into account the fact that the complete free 
volume can not be filled with drops. Note that, for different symmetries of the suspended particles, 
the value of the constant c will be different. This is specially important in the case of drops which 
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may change their shape when subjected to strong shears. Therefore, c contains information about 
the maximum packing of drops the system may allocate. 

Therefore, if excluded volume effects are taken into account, the suspended phase contribution 
to the stress tensor is given by 

iff re. _ N r ^{l),y 

^^'^ -y- cNVd J ^^'^ 

For finite-sized droplets, this relation leads to the result that the Taylor expression scales with the 
excluded volume factor — ccf)) instead of (p, and thus gives the following expression for the 
viscosity of an emulsion 

1 + 2.5K 



Vr (0) 



1 + 



(8) 



l + K \l-c 

Taylor's formula is then recovered in the limit of very small volume fractions {(p — > 0). Eq. ^ 
explains the experimental fact that the viscosity diverges at volume fractions </><!, implying that 
this effect is related to the excluded volume effects. Similar arguments used in the case of solid 
particles lead to excellent agreement between experiments and theory, see Ref. 1^. 

The central result of the present analysis is Eq. ([8]) , which shows that the correct scaling factor 
describing the dependence of the viscosity of an emulsion of finite-sized droplets on the filling 
fraction is the effective filling fraction 

= T^y 

where c is a structural or crowding factor that takes into account the arrangements of the droplets 
in the emulsion and is given by 

0='—^. (10) 

<Pc 

where (pc is the filling fraction at the divergence. The importance of (peff is that it incorporates 
correlations between particles by taking into account geometrical information having the charac- 
teristic that 4>eff ^ (p &t low volume fractions and (peff = 1 at maximum packing. This quantity 
plays an important role in determining the correct dependence of the viscosity of the emulsion, as 
we will show next. 

At low volume fractions, the shear viscosity r] of an emulsion is given by Eq.dSj). In order 
to extend it to larger filling fractions (p, we will use a differential effective medium approach in 
which a concentrated emulsion is obtained from an initial continuous phase by successively adding 
infinitesimally small quantities of droplets to the system until the final volume fraction of the 
dispersed phase is reached. 
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The usual implementation of the differential theory to obtain a concentrated emulsion up to a 
given stage, consists in treat it as a homogeneous effective medium of viscosity r] {(j)) into which we 
add a quantity Ai;^ of new droplets. Then, the new viscosity ri[(j) + Sep) is calculated using Taylor's 
expression to give 

,(^ + i^) = „w(l + ''<ti^f"-' A^). (11) 

Here, it is important to notice that in order to allocate the new droplets maintaining constant 
the volume of the system, one has to remove part of the effective medium (1 — A(^) in order 
to allocate the new quantity (A0) of droplets. Therefore, the new filling fraction is given by 
+ 50 = 0(1 — A(/>) + A0, from which one finds 

A,^J^. ,12) 

A shortcoming of using Eq. (jlip and the bare filling fraction is that it does not contains the 
correlations generated by the excluded volume effects, that is, it assumes that all the volume of 
the emulsion before new droplets are added is available to the new droplets. However, as we have 
mentioned previously, this is not possible since the droplets can occupy only the free volume and 
therefore scale according to 0e// = 0(1 — c^)""*^, Eq. Q. 

These considerations suggest that the differential theory must be performed in terms of the 
scaling variable 0e// and not on the bare filling fraction itself. 

Therefore, making the substitution of by 0e// in Eqs. ijllh and (I12p and integrating we obtain 

, ^ 2r?,(0e//) + 5j^ V/^ 5/2 
^'^"■^■^^ V 2T5^ ) " ^ "^^^^^ ' ^ ' 

that, in terms of the bare filling fraction takes the form 

Expressions ()13p and (I14p are similar to other ones based on the differential method, see, for 
example, Eqs. (l2])-((4l). However, our model crucially makes use of the geometrical information of 
the system from the start through 0e//. This makes an essential difference since, in contrast to the 
usual differential procedure, the use of 0e// as the integration variable incorporates correlations 
between droplets of the same recursive stage, resulting in an improvement of the quantitative 
description when compared to experimental data, as will be shown in the next section. 
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III. RESULTS 

In Figure 1 we show the relative viscosities predicted form the improved DEMT, Eq. ()13p . and 
Pal's model 2, Eq.dS]), as a function of the viscosity ratio K. Both models initially show a nearly 
constant value followed by an increase in the range 10^^ < K < W^, approximately, and finally 
becomes nearly constant again for large values of K. They also predict an increase with increasing 
(j). In this figure, 0c was taken to be 0.7404. A comparison between both models show that the 
improved DEMT, Eq. (jl3p . predicts higher values of the relative viscosity than Pal's model 2, Eq.([3]), 
over the full range of K. The difference between both models is larger for larger 0. 

The various data sets considered for comparison with the models are the same as the ones used 



in Ref. [8(]. They consist on emulsions covering a viscosity ratio 3.87 x lO""^ < K < 3.25 x 10^ 
and with the characteristic that the capillary number is small and that the emulsions are stable 
(unfloculated) . 

Figure 2 shows comparisons between the experimental data and predictions of improved DEMT, 
Eq. (jl3p . and Pal's model 2, Eq.Q, for different values of the viscosity ratio K. As can be seen, 
improved DEMT describes the experimental data slightly better than Pal's model 2. However, 
the more significant difference between both models is the fitting parameter (j)c which varies for 
different values of K. In all cases the improved DEMT predicts larger values of (pc- Figure 3 
shows the independence of the relative viscosity on the mean droplet size. The comparison with 
the theoretical models show again a better fit for the improved DEMT. In Fig. 4 we compare 
experimental relative viscosity data at a fixed (j) = 0.5 as a function of the viscosity ratio K. Both 
models show a good agreement with experimental data, nonetheless approximate fitting values of 
(f)c are very different in both models, being larger for the improved DEMT. 

According to Eq.([T3]), the viscosity, plotted in the form ^r {(p)'^^^ [(2??r (0) + 5K)/{2 + 5K)]~^/^ 
versus should be system independent. This is confirmed in Fig. 4 which shows all the viscosity 
data, and compare them with the improved DEMT prediction. Pal's model 2 can not be plotted 
in this sample independent form since it produces different curves for different values of 0c. As can 
be seen, the agreement between data points and improved DEMT is excellent. 

IV. CONCLUSIONS 

An improved viscosity equation is developed for emulsions of spherical droplets starting from 
the Taylor's equation for the viscosity of very dilute emulsions HI] and using a differential effective 
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medium approach. Geometrical information about the packing of the droplets is included through 
an effective volume fraction, (j)eff, which approaches (p at low concentrations and becomes one at 
the critical concentration. The model improves the results obtained with Pal's model 2, Ref. ^ 
for the relative viscosity, describes experimental data very well, and reduces to Taylor's expression 
at low concentrations. In the limit X ^ oo it reduces to the expression recently found for a 
suspension of hard spheres . 
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Figure Captions 

Fig. 1. Relative viscosity rjr ((f)) = rj {(f)) /r/o predicted from the improved DEMT, Eq. (jl3p (solid 
line) and Pal's model 2, Eq.Q (dashed line), (pc is taken to be 0.7404. 

Fig. 2. Comparisons between the experimental data [Pal (2001)] and best fits of the theoretical 
predictions for the volume dependence of the relative viscosity for different values of the viscosity 
ratio K. The solid line is the result of the improved DEMT and the dashed line the prediction of 
Pal's model 2. 

Fig. 3. Comparisons between the experimental data and predictions of improved DEMT (solid 
line), Eq. (jl3p and Pal's model 2 (dashed line), Eq.Q. The best fit for the improved DEMT was 
obtained with (pc = 0.67 and for Pal's model 2 0c = 0.61. 

Fig. 4. Viscosity ratio (K) dependence of the relative viscosity. The solid line is the result of 
the improved DEMT with (f>c = 0.7404 and the dashed line is the prediction of Pal's model 2 with 
(f)c = 0.63. 

Fig. 5. Master curve of 1 — 4>eff consisting of all the experimental viscosity data. The solid line 
is the result of the improved DEMT. 
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